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ABSTRACT 

We investigate the impact of photochemistry and X-ray ionization on the molecular composition 
of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a 
sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and 
X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical 
network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recal- 
culate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent 
reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental compo- 
sition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular 
composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, 
HCO+, N 2 H+, H 2 0, C0 2 and CH 3 OH. The only molecule significantly affected by the X-ray ion- 
ization is N 2 H + indicating it is safe to adopt existing approximations of the X-ray ionization rate in 
typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of 
neutral molecules in the outer disk, highlighting the importance of taking into account the shape of 
the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the 
gas-phase chemistry due to the adjustment of the H/H 2 and C + /C ratios. The disk ionization frac- 
tion is not significantly affected by the methods adopted to calculate the photochemistry and X-ray 
ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a 
layer, Z/R < 0.1 - 0.2, in the disk midplane, within R w 200 AU. 

Subject headings: astrochemistry - protoplanetary disks - stars:formation - ISM:molecules 



1. INTRODUCTION 

Protoplanetary disks have several vital functions in 
star and planet formation: they (I) aid the dissipation 
of angular momentum away from the protostellar sys- 
tem, (2) allow the efficient accretion of matter from the 
constituent cloud material onto the central star, and (3) 
contain the material which may eventually form an ac- 
companying planetary system. 

Protoplanetary disks are ph ysically and thus, chemi- 
cally, complex objects (see e.g. jBergin et~ani2007h . They 
are heavily irradiated by UV radiation from their par- 
ent T Tauri star and are permeated by X-rays and ex- 
cess UV photons thought to arise from an accretion 
shock generate d as disk material impinges upon the 
stella r surface (|Herbig fe Goodrich] 119861 : iKastner et al.l 
119971) . Beyond a radius r > 100 AU, the UV radi- 
ation field originating from the parent star decreases 
in strength due to a combination of absorption of UV 
photons by the intervening disk material and geomet- 
rical dilution. Here, irradiation by UV photons orig- 
inating from the interstellar radiation field (ISRF) in- 
creases in importance and, as a result, the wavelength 
dependence (or shape) of the radiati on field varies as 
a function of disk radius and height (|Aikawa fc Herbst] 
119991 I Willacv fc Langerl l200dh . As a result, in mod- 

catherine.walsh@qub.ac.uk 

1 Astrophysics Research Centre, School of Mathematics and 
Physics, Queen's University Belfast, University Road, Belfast, 
Northern Ireland, UK, BT7 INN 

2 Department of Astronomy, Graduate School of Science, Ky- 
oto University, Kyoto 606-8502, Japan 

3 Department of Earth and Planetary Sciences, Kobe Univer- 
sity, 1-1 Rokkodai-cho, Nada, Kobe 657-8501, Japan 



ern disk models, wavelength-dependent radi ative trans 
fer is preferred (|van Zadelhoff et al.l I2003D as i s the 
inclusion of UV excess radiation dBergin et al.l 1200 
Nomura fc Millar! 12005ft and X-ravs (|Glassgold et al 



19971: lAikawa fe Herbstlll999| ). Along with cosmic-rays 



the X-ray and UV radiation controls the ionization frac- 
tion in the disk which has consequences on the disk 
accretion rate and the location and extent of 'dead 
zones', regions where accretion is potent ially suppressed 
(jBalbus fc Hawlevl IToMl : lOammiel fl996h . The varying 
radiation field will also have a direct effect on the disk 
chemical structure, controlling the abundance and dis- 
tribution of atoms, ions and molecules through pho- 
tochemistry and influencing the molecular composition 
of the icy grain mantle via non-thermal desorption 
and this has been demonstrated in many works (see 
e.g.. Ivan Zadelhoff et ail 120031: lAikawa fc Nomural 120061: 
Wals h et al.l 120101: iKamp et al.l 120101: iVasvunin et al.l 



20111) . For these reasons, the treatment of photo pro 



cesses in protoplanetary disks should be thoroughly in- 
vestigated, in order to aid the interpretation of obser- 
vational data, especially with the impending completion 
of the Atacama Large Millimeter Array (ALMA) which, 
for the first time, will enable the observation of molecu- 
lar emission from nearby (~ 140 pc) protoplanetary disks 
on around sub-milli-arcsecond scales with unprecedented 
spectral resolution. 

A plethora of molecules have been detected in proto- 
planetary disks via line emission in the (sub)mm and in- 
frared regions of the electromagnetic spectrum. Early ob- 
servations at (sub)mm wavelengths were made using the 
James Clerk Maxwell Telescope ( JCMT) and IRAM 30 m 
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radio telescope (e.g.. iKastner et al.l |1997 | : Dutrev et al.l 



[19971 Ivan Zadelhoff et al.1 120011: iThi et al.lT2004D . with 
more recent detections using the Submillimeter Array 
(SMA) (e.g., IQi et alJ I20M [20081: IQberg et al.1 l2010h . 
Most molecules observed in this spectral region are small, 
simple, abundant molecules, molecular ions and radicals 
(e.g., CO, CN, CS, HCO+, N 2 H+, HCN) and associated 
isotopologues (e.g., 13 CO, DCO+, and C 34 S). A recent 
survey of disks around T Tauri and Herbig Ae stars has 
led to the fi rst successful detec tion of SO in a circum- 
stellar disk (jFuente et al.1 12010). with the authors also 
reporting a tentative detection of H 2 S. Due to the limi- 
tations of existing telescopes and the small angular size 
of disks, the most com plex molecule observed to date is 
form aldehyde, H 2 CO (|Dutrev et all 119971; lAikawa et all 
l200l . 

Use of the Infrared Spectrograph (IRS) on the 
Spitzer Space Telescope increased the inventory of 
gas-phase molecules detected in disks to include OH 
H 2 Q CQ 2 and CyH? (lLahuis et al.ll2006HCarr fc Naiital 
120081: iSalvk et al.1 120081: iPontoppidan et all I2010D . Ex- 
isting (sub)mm observations probe the colder, outer 
disk material whereas infrared observations probe 
the warmer gas in the inner disk surface in the 
so-called 'planet-forming' region. There have also 
been detections of water ice absorption features in 
'edge-on' T Tau r i systems (iCreech-Eakmanl 120021 : 
iTerada et all [20071: iSchegerer fe Wolj 120101). More re- 
centlv. iHogerheijde et al.l 1 20111 ) report the first detec- 
tion of the ground-state rotational emission lines of both 
spin isomer states of water in a protoplanetary disk using 
the Heterodyne Instrument for the Far-Infrared (HIFI) 
mounted on the Herschel Space Observatory. These sets 
of observations, of both gas and ice, give us a reason- 
ably sufficient benchmark with which we can compare 
our results. 

In iWalsh et al.l ([2010D , henceforth referred to as Pa- 
per I, we used t h e phy sical disk model de s cribed in 
INomura fc Millarl ([20051 ) and iNomura et all ([2007[ ) to 
compute the chemical structure of a typical protoplan- 
etary disk on small scales (sub-milli-arcsecond in the 
inner disk for an object at the distance of Taurus, 
~ 140 pc), investigating the effects of the addition of 
non-thermal desorption mechanisms (cosmic-ray-induced 
desorption, photodesorption and X-ray desorption) and 
grain-surface chemistry on the disk chemical structure. 
In that work, we presented results from models in which 
we approximated the photoreaction rates by scaling the 
rates f rom the UMIST Dat abase for Astrochemistry or 
UDf A (jWoodall et al.ll2007|) , which assume the interstel- 
lar radiation field (ISRF), by the wavelength integrated 
UV flux at each point in the disk (see Section 12.2.11 for 
further details). Here, we report results from models in 
which we explicitly calculate the photodissociation and 
photoionization rates taking into consideration the UV 
spectrum at each point and the wavelength-dependent 
absorption cross section for each photoreaction. In addi- 
tion, we recalculate the X-ray ionization rate everywhere 
in the disk accounting for the elemental composition of 
the gas and include the direct X-ray ionization of ele- 
ments, in both cases, using the X-ray energy spectrum 
at each point. 

In Section 12.11 we give a brief overview of our physi- 



cal model. In Section 12721 we describe our chemical net- 
work and processes we include in our calculation of the 
chemical structure with a thorough description of the 
methods used to compute the photochemical and X-ray 
ionization rates (Sections 12.2.11 and I2.2.2) respectively). 
In Section l2~3l we describe the theory behind the identifi- 
cation of regions of our disk in which angular momentum 
transport and thus, accretion, may be suppressed. The 
results of our calculations are presented in Section[3]with 
a summary given in Section [4l 

2. PROTOPLANETARY DISK MODEL 
2.1. Physical Model 
Th e physical model of a pr otoplanetary disk we use is 



from INomura fc Millarl ([2005D with the addi tion of X-ray 
heating as described in lNomura et al.1 (|2007l ) . The degree 
of ionization in the disk depends on the disk parameters 
adopted and the resulting surface density distribution. 
The theoretical foundation of ou r model comes from the 
standa rd ac cretion d isk m odel of iLvnden-Bell fc Pringld 
(|1974l ) and iPringld ([19811 ) which defines a surface den- 
sity distribution for the disk given the central star's 
mass and radius and a disk accretion rate, M. The 
kinematic viscosity in the disk is param e terise d accord- 
ing to the work of iShakura fc Sunvaevl ([1973f L the so- 
called, a-prescription. We consider an axisymmetric 
disk surrounding a typical T Tauri star with mass, 
M* = 0.5 M , radius, i?* = 2 i? and effective tem- 
perature, T* = 4000 K. We adopt a viscous parameter, 
a = 0.01 and a mass accretion rate, M = 10 -8 M yr . 

We use a model X-ray spectrum created by fitting the 
observed XMM-New ton spectrum o f the cl assical T Tauri 
star, TW Hydrae (jKastner et al.l |2002) with a two- 
tempera ture thin thermal pla sma model (MEKAL model 
see e.g., iLiedahl et al.l (|1995D ). The derived best- fit pa- 
rameters for the plasma temperatures are kT\ = 0.8 keV 
and fcT 2 = 0.2 keV and for the foreground interstellar 
hydrogen column density, Ah = 2.7 x 10 20 cm~ 2 . For 
the X-ray extinction, we include attenuation due to ion- 
ization of all elements and Compton scattering by hy- 
drogen. The X-ray luminosity is Lx ~ 10 30 erg s _1 and 
the resulti ng high-resolution X-r ay spectrum is given in 
Figure 1 of INomura et alJ (|2007l ) assuming a distance to 
source of 56 pc, and is reproduced in binned form here 
in the right-hand panel of Figure [TJ 

The UV radiation field in disks has two sources, the 
star and the interstellar medium. In this disk model, 
the radiation field due to the T Tauri star has three 
components: black-body emission at the star's effective 
temperature, optically thin hydrogenic bremsstrahlung 
emission and strong Lyman-a line emission. For the 
UV extinction, we include absorption and scattering by 
dust grains. We assume the dust and gas in the disk 
is well mixed and adopt a dust-size distribution model 
which reprod uces the observational extinc tion curve of 
dense clouds (fWcingar tner fc Dra inc 200fl). The calcu- 
lation of the dust opacity in the disk is described in 
Appendix D of INomura fc Millarl (|2005[ ) with the result- 
ing wavelength-dependent absorption coefficient shown 
in Figure D.l. The total FUV luminosity in our model 
is Lyv ~ 10 31 erg s _1 with the calculation of the radia- 
tion field in the di s k desc ribed in detail in Appendix C of 
INomura fc Millarl (|2005| ). We display the resulting stel- 
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lar flux density in the disk surface at a radius of 1 AU, 
including each individual component, in the left-hand 
panel of Figure [TJ The main source of UV photons short- 
ward of 2000 A is due to Brehmsstrahlung and Lyman-a 
radiation with the Lyman-a line contributing around 10 3 
times the UV continuum photon flux at ks 1216 A over an 
assumed FWHM of « 2 A (see e.g. JHerczeg et al.ll2002T ). 

The resulting disk physical structure is given in Fig- 
ures 1 and 13 in Paper I and we refer readers to the 
Appendix of that publication for a thorough discussion. 
In Figure [21 we display the gas and dust temperatures in 
K (top right), the gas number density in cm -3 (top left), 
the wavelength-integrated UV flux (bottom left) and X- 
ray flux (bottom right) both in units of erg cm -2 s , 
as a function of disk radius and height (scaled by the ra- 
dius). In the temperature plot, the colour map represents 
the gas temperature whereas the contours represent the 
dust temperature. As expected, the disk surface closest 
to the parent star is subjected to the largest flux of both 
UV and X-ray radiation. The disk midplane is effectively 
completely shielded from UV photons over the radial ex- 
tent of our disk model. The higher energy X-ray photons, 
although resulting in a lower flux in the disk surface, are 
able to penetrate the disk more effectively leading to a 
small, yet appreciable, X-ray flux in the disk midplane 
beyond w 10 AU. 

2.2. Chemical Model 

As in Paper I, our gas-phase chemistry is ex- 
tracted from the latest release of the UMIST 
Database for Astroche mistry or 'Rate06', a vailable at 
|http : //www . udf a . net] (|Woodall et al.l l2007t) . We use 
almost the entire Rate06 network removing only those 
reactions involving fluorine- and phosphorus-containing 
species. Here, we also include the subset of three- 
body reactions from Rate06, since these may play a sig- 
nificant role in the densest regions of our disk model 
(n ~ 10 15 cm -3 ). 

We allow gas-grain interactions i.e., the accretion of 
gas-phase species onto dust grains with removal of the 
grain mantl e via thermal desorption, cosmic-ray-induced 
desorption (|Leger et al.lll985HHasegawa &: Herbstill993D 
and p hotodesorptio n (|Westlev et al.l 119951 : lOberg et al.l 
120071: lWillacvll2007l ). In Paper I we provided a thor- 
ough description of the methods used to determine our 
accretion and desorption r ates. We also i nclud e the 
grain-surface network from iHasegawa et al.1 (|1992f ) and 
lHasegawa fc HerbstJ (| 19931) . 

Here, we calculate our photodesorption rates taking 
into account the molecular composition of the grain man- 
tle. Experiments into the photodesorption of UV ir- 
radiated ices conducted by lOberg et al.l (2009a. bj) sug- 
gest that the photodesorption yields are dependent on 
the ice composition with pure CO, Na, CO2 and H2O 
ices giving different desorption yields (listed in Tabled]). 
The picture is further complicated by evidence of codes- 
orption in mixed ices, with the photodesorption yield of 
N2 increasing when present in a 1:1 N2/CO ice mixture 
(jOberg et al. 2007). More recent work suggests that in 
organic ices, irradiation by UV photons initiates chem- 
istry with photodesorption products detected other than 
those present in the original ice mixture (jOberg et al.l 
l2009d ). These latter effects are not included in this model 



TABLE 1 
Photodesorption Yields 



Species 


Yield (molecules photon 1 ) 


Reference 


CO 


2.7 x icr 3 


1 


N 2 


1.8 x icr 4 


1 


C0 2 


2.3 x icr 3 


2 


H 2 


1.3 x icr 3 


2 


All other species 


3.0 x 10~ 3 


3 



Refe ren ces. — (1) lOberg et all I|2009al 1. (2) lOberg eTal 
q2009br i. (3) IWestlev et alJ l|1995l l 



TABLE 2 
Molecular Binding Energies 



Species 


Binding Energy (K) 


Reference 


CO 


855 


1 


N 2 


790 


1 


CO2 


2990 


2 


C 2 S 


5320 


2 


H 2 


4800 


3 


SO2 


5330 


2 


NH 3 


2790 


3 


CH 4 


1090 


1 


HCOOH 


5000 


5 


CH3OH 


4930 


3 


CH3CHO 


3800 


5 


C 2 H 6 


2300 


5 


HCOOCH3 


4000 


5 


CH3OCH3 


3300 


5 


C 2 H 5 OH 


5200 


5 


References 


— (1) lOberg et al.l 


(2005), (2) 



IHerrero et"aLl (t20Tol 1 . (5) lOberg et aT] d2009d 1 

for simplicity, however, they should be investigated in fu- 
ture models. 

In Paper I, we used a constant photodesorption yield 
of 3 x 10~ 3 phot o n -1 a s determined for pure water ice 
by Wcstlc v" et al.l ((1995). To take into consideration the 
composition of the grain mantle we now determine the 
photodesorption rate for a specific species according to 



fl 



where Guv is the radiation field in units of photons 
cm~ 2 s _1 , Y^y photon -1 is the specific desorption yield 
of species i, as listed in Table Q] x\ — nl/nl ot is the 
fractional abundance of species i on the dust grains and 
o~a and Xd are the dust-grain geometrical cross section 
and fractional abundance, respectively. We have also re- 
viewed our set of molecular desorption energies in light of 
recent experimental results and in Table [5] we list those 
species for which the binding energies, E^, have been 
updated. 

Our initial fractional abundances are the result of 
a dark cloud model run with typical molecular cloud 
parameters i.e., T = 10 K, n(Ha) = 10 4 cm -3 and 
A v — 10 mag using the set of o xygen-rich low - metall icity 
elemental abundan ces from iGraec el et al.l (|1982| ) as 



2007). In the gen- 



listed in Table 8 of iWoodall et al.1 
eration of our initial abundances we allow for freeze out 
and thermal desorption and we extract abundances at a 
time of 10 5 years which is thought to be the age of dark 
clouds on the brink of star formation. 
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Fig. 1. — Stellar UV flux density in the disk surface at 1 AU in photons cm" 
cm -2 s — 1 keV — 1 (right). The latter assumes a distance to source of 56 pc. 
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2.2.1. Photochemistry 

The photoreaction rates in Rate06 are calculated as- 
sumin g the UV ra diation field is given by the Draine 
field (jDraind [1978), an adequate assumption for the 
unshielded ISRF. For use in chemical models, these 
rates are subsequently parameterised according to optical 
depth or A v . In Paper I, we approximated our photore- 
action rates in the disk, /c p h, by scaling the rates from 
Rate06 using the wavelength-integrated UV flux at each 



point, Guv (7, z) 



r 2000 a „ , x , . . 
hi2°A Gx{r,z)d\i.e., 



G 



UV 



Go 



k s" 



(2) 



where Go is the unshielded interstellar UV flux and 
ko is the rate calculated for the unshielded interstellar 
medium. Note, G\(r, z) includes both the stellar and in- 
terstellar components of the radiation field. This approx- 
imation is unsuitable for use in protoplanetary disks as 
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these objects are irradiated by UV radiation with three 
components: the central star, UV excess due to accre- 
tion and the interstellar radiation field. Hence, the UV 
spectrum at each point in the disk will not only vary 
with disk radius and height, but will also bear no resem- 
blance to the interstella r radiation field (see Figure 4 in 
iNomura fc Mi llar (|2005ft ). For these reasons, we have in- 
vestigated a recalculation of the photorates in the disk 
taking into consideration the UV spectrum at each point 
van DishoeckllT987l : Ivan Dishoeck et al.ll2006ft . 

The photodissociation rate due to continuous absorp- 
tion, fcp h , is calculated using 



h c — 



a(A)J(A) dA s" 



(3) 



where, A is the wavelength, tr(A) is the cross section and 
/(A) is the mean intensity of UV radiation. The rate for 
photoionization, fc p i, is calculated using the same equa- 
tion. For indirect photodissociation via absorption into 
a bound upper state, u, from a lower state, I, the rate is 



(4) 



where, A u i is the wavelength of the line transition, / u i 
is the oscillator strength of the transition and ^t u i is 
an efficiency factor. The parameters, e, c and m are 
the electron charge, the speed of light and the atomic 
or molecular mass, respectively. The total photodis- 
sociation rate is found by summing over all possible 
channels. The pho toreaction cross sect io ns, tr (A), are 
those adopted in Ivan Zadelhoff etafl (12003ft w hich 
originate from calculations by Ivan Dishoeck (11 9881). up- 
dated by Uansen et al.l (|1995al fbl) and Ivan Dishoeck et al.l 
(12006ft. The photo cross sections are dow nloadable from 
http : / / www . strw . leidenuni v.nl/[ ^ewine/pho to/[ 
For species which do not have a calculated cross section, 
we use the rate for a similar type of molecule. 

We should note here we include Lyman-a radiation 
in the calculation of the wavelength-integrated UV flux 
(Equation [2]) which ac counts for approxim ately 85% of 
the total flux (see e.g.. iBergin et all 12003ft . however, we 
ignore it in the calculation of the wavelength-dependent 
UV radiation field and subsequent photorates (see Equa- 
tions [3] and 3]). The exclusion of Lyman-a is primar- 
ily due to the difficulties in treating the scattering of 
Lyman-a photons from the surface into the disk (see e.g., 
IBergin et al.ll2003ft . Lyman-a scattering differs from that 
of the background UV photons since the scattering oc- 
curs predominantly by H atoms rather than dust grains. 
Indeed, Lyman-a radiation has been historically ne- 
glected in protoplanetary disk models and has only very 
recent ly been addressed in the work bv lBethell &: Berginl 
(|2011ft . However, the photorates calculated according 
to Equations [3] and 2] are more accurate than those de- 
termined using Equation [5] since the shape of the back- 
ground radiation field is included in the calculation. The 
photorates calculated using Equation [5] include the total 
flux of UV photons (background plus Lyman-a) heav- 
ily overestimating the strength of the UV field at wave- 
lengths other than the Lyman-a wavelength (« 1216 A). 
The effect of including the wavelength-dependence in the 
photorates is apparent in differences in the relative abun- 
dances of molecules (see Table|3j) e.g., the column density 



ratio of CO2/H2O changes from w 0.045 in model UV- 
old to ~ 1.90 in model UV-new at 100 AU. We discuss 
the possible effects of including Lyman-a radiation in the 
calculation of the photorates in Section [3761 

We include the self- shielding of Ho using the approx- 
imation from iFederman et al.1 ([1979) in the generation 



of our physical model giving us the initial conditions 
in our disk, however, we do not explicitly include the 
self- and mutual shielding of H 2 and CO in our calcula- 
tion of the subsequent chemical structure. We find the 
dominant component of the radiation field in the disk 
surface is the radial component which is the direct stel- 
lar radiation (the vertical component consists of both 
the contribution from the interstellar radiation field and 
scattered stellar radiation). In addition, throughout the 
majority of the disk, the stellar radiation dominates over 
the interstellar radiation. Hence, we argue against the 
validity of adopting the usual plane-parallel approxima- 
tion for the calculation of the self-shielding factors and 
the application of shielding factors computed for irra- 
diated interstellar clouds, to protoplanetary disks. To 
correctly include the effects of self- and mutual shielding 
in disks, a self-consistent two-dimensional treatment is 
needed which takes into consideration the time-varying 
H 2 and CO abundances throughout the disk, the col- 
umn densities in the radial and vertical direction and 
the two-dimensional physical structure of the disk which 
will effect the line widths and line strengths and hence, 
shielding factors. We discuss this issue further in Sec- 
tion [3751 



2.2.2. X-ray Ionization 

The model of INomura et al.l ([2007ft calculates an over- 
all X-ray ionization rate at each point in th e disk accord- 
ing to the theory of iMalonev et al.l (|1996ft . The rate at 
each point is calculated assuming a power-law fit for the 
X-ray absorption cross section dependent on X-ray en- 
ergy and it is these approximate rates that are used in 
our chemical calculations in Paper I. In this work, we re- 
calculate the X-ray ionization rate everywhere in the disk 
taking into account the X-ray energy spectrum, Fx{E), 
at each point and the explicit elemental composition of 
the gas (|Glassgold et al.lll997ft . The overall X -ray ion- 
ization rate, Cxr, is given by summing over all elements, 



c 



XR 



E 



x k a k (E)F x (E) 



E-E k 
Ae 



dE: 



( 5 ) 

where, for each element, k, E k is the ionization potential, 
<r k {E) is the cross section and x k is the fractional abun- 
dance with respect to H nuclei density. In this expression, 
the number of secondary ionizations per unit energy pro- 
duced by primary photoelectrons, N scc , is given by the 
expression, (E — E k )/Ae, where Ae is the mean energy 
required to make an ion pair (« 37 eV). Only the number 
of secondary ionizations needs to be considered as this is 
generally much larger than the number of primary ion- 
ization events. Typically, each keV of secondary electron 
energy produces an average of 1000/37 « 27 ion pairs so 
that N scc 3> N pi i. Note that X-rays interact only with 
atoms, regardles s of whether an atom is bound within a 
molecule or free ([Glassgold et al.|[l997ft . 
We have also added the direct X-ray ionization of ele- 
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TABLE 3 
Column Densities 



Species UV-old XR+UV-old UV-new XR+UV-new 



0.1 AU 



H 


1.4(21) 


1.4(21) 


4.5(20) 


4.2(20) 


H 2 


1.1(26) 


1.1(26) 


1.1(26) 


1.1(26) 


CO 


9.7(21) 


9.5(21) 


9.5(21) 


9.6(21) 


HCO+ 


2.1(14) 


2.5(14) 


3.2(12) 


1.7(12) 


HCN 


5.2(19) 


5.0(19) 


5.4(19) 


4.9(19) 


CN 


8.2(13) 


6.9(13) 


4.3(13) 


6.5(13) 


cs 


1.3(17) 


1.5(17) 


1.5(17) 


1.5(17) 


C 2 H 


1.0(15) 


9.9(14) 


9.6(14) 


9.8(14) 


H 2 CO 


6.9(15) 


8.2(15) 


6.5(15) 


8.3(15) 


N 2 H+ 


3.7(09) 


8.7(08) 


2.0(11) 


9.2(10) 


OH 


6.0(16) 


6.0(16) 


6.1(15) 


5.1(15) 


H 2 


2.8(22) 


2.8(22) 


2.8(22) 


2.8(22) 


co 2 


3.8(18) 


5.4(18) 


3.7(18) 


5.5(18) 


C 2 H 2 


6.1(18) 


9.6(18) 


9.3(18) 


9.1(18) 


CH 3 OH 


2.3(18) 


2.4(18) 


2.3(18) 


2.4(18) 



1 AU 



H 


1.1(21) 


1.1(21) 


8.8(20) 


8.3(20) 


H 2 


1.9(25) 


1.9(25) 


1.9(25) 


1.9(25) 


CO 


1.9(21) 


1.9(21) 


1.9(21) 


1.9(21) 


HCO+ 


1.7(14) 


1.9(14) 


9.4(13) 


6.5(13) 


HCN 


2.0(18) 


2.0(18) 


2.0(18) 


2.0(18) 


CN 


2.7(14) 


2.7(14) 


1.0(14) 


1.3(14) 


CS 


8.7(12) 


8.4(12) 


6.2(12) 


6.3(12) 


C 2 H 


4.1(14) 


4.4(14) 


9.1(13) 


1.5(14) 


H 2 CO 


3.1(13) 


3.1(13) 


3.2(13) 


3.1(13) 


N 2 H+ 


5.2(10) 


9.3(09) 


2.7(12) 


5.1(11) 


OH 


1.7(16) 


1.8(16) 


2.4(16) 


2.2(16) 


H 2 


2.2(21) 


2.2(21) 


2.2(21) 


2.2(21) 


co 2 


7.7(20) 


7.7(20) 


7.7(20) 


7.7(20) 


C 2 H 2 


1.6(14) 


1.9(14) 


1.5(14) 


1.8(14) 


CH3OH 


1.0(15) 


1.0(15) 


1.0(15) 


1.0(15) 



10 AU 



H 


7.3(20 


7.2(20) 


6.1(20) 


5.7(20) 


H 2 


2.6(24) 


2.6(24) 


2.6(24) 


2.6(24) 


CO 


2.0(20) 


2.0(20) 


2.0(20) 


2.0(20) 


HCO+ 


4.8(13) 


4.5(13) 


1.4(14) 


1.6(14) 


HCN 


8.0(14) 


2.1(15) 


4.5(14) 


6.9(14) 


CN 


7.3(13) 


8.1(13) 


3.0(14) 


3.4(14) 


CS 


4.7(13) 


4.4(13) 


4.3(13) 


5.5(13) 


C 2 H 


5.4(13) 


6.7(13) 


1.1(14) 


1.6(14) 


H 2 CO 


2.0(12) 


2.0(12) 


1.2(12) 


8.1(11) 


N 2 H+ 


2.0(10) 


5.7(09) 


2.8(11) 


9.2(10) 


OH 


8.9(15) 


9.4(15) 


2.1(16) 


2.3(16) 


H 2 


5.4(15) 


5.9(15) 


8.5(16) 


9.2(16) 


C0 2 


2.8(16) 


3.7(16) 


2.7(17) 


6.5(17) 


C 2 H 2 


7.5(13) 


3.2(13) 


7.9(13) 


3.7(13) 


CH3OH 


1.0(08) 


9.1(07) 


5.3(09) 


5.2(09) 



100 AU 



H 


3.3(19) 


1.7(19) 


3.2(19) 


1.5(19) 


H 2 


2.0(23) 


2.0(23) 


2.0(23) 


2.0(23) 


CO 


9.1(18) 


9.7(18) 


8.9(18) 


9.3(18) 


HCO+ 


3.9(13) 


3.6(13) 


5.8(13) 


4.3(13) 


HCN 


4.8(14) 


2.5(14) 


2.3(14) 


2.1(14) 


CN 


2.4(14) 


2.6(14) 


5.0(14) 


4.1(14) 


CS 


3.1(13) 


3.7(13) 


1.8(13) 


1.7(13) 


C 2 H 


1.0(14) 


1.2(14) 


7.7(13) 


5.4(13) 


H 2 CO 


4.4(12) 


3.7(12) 


8.2(12) 


2.9(12) 


N 2 H+ 


1.9(11) 


8.8(10) 


1.5(12) 


5.5(11) 


OH 


9.0(14) 


3.0(14) 


6.3(15) 


3.3(15) 


H 2 


4.9(16) 


2.5(16) 


3.1(16) 


1.4(16) 


C0 2 


2.2(15) 


2.1(15) 


5.9(16) 


8.8(16) 


C 2 H 2 


3.0(13) 


3.1(13) 


2.4(13) 


1.6(13) 


CH3OH 


4.1(10) 


4.0(11) 


7.9(10) 


2.1(11) 


Note. - 


a (ft) means 


a X 10 b 







ments, the rate for which is given by 

Cfc = / a k (E)F x (E) dEs- 1 . (6) 

usin g the ionization cro ss sections for each element, k, 
from lVerner et alj (|1993l ). 

2.3. Disk Ionization Fraction 

In addition to investigating the importance of photo- 
chemistry and X-ray chemistry in protoplanetary disks, 
we have determined the location and extent of potential 
dead zones where accretion may be inhibited. Angular 
momentum transport in disks is thought to arise from 
turbulence initiated by a weak-field magne torotational 
instability or MRI (Ba lbus fe Hawlevlll99lD . hence, ac- 
cretion may be inefficient in regions where the instability 
is suppressed. The turbulence generated by the insta- 
bility can sustain a disordered magnetic field to which 
the gas is coupled. The degree of the coupling, in turn, 
depends on the ionization fr action in the disk. 

Following [Gammic (1996), we can define a magnetic 
Reynolds number, ReM, everywhere, 



where Va ~ a 1 ^ 2 c s is the Alfven speed, a function of a, 
the scaling parameter for th e viscosity (v = ac.^H) from 
the accretion disk model of iShakura fc Sunyaevl (|1973f l 
and c s , the sound speed in the disk. H is the disk scale 
height given by H = c s /0, where f2 is the Keplerian 
velocity at a particular radius, R, and rj is the magnetic 
resistivity which is related to the electron fraction, x = 
n c /nn by 

n = (6.5 x lO 3 )*" 1 cmV 1 , (8) 

(see iGammiel (|1996h for further details). Accretion is 
likely suppressed in regions where the magnetic Reynolds 
number, ReM, falls below a critical value, i?e^J lt . This 
parameter determines the degree to which the ionized 
gas is effectively coupled to the magnetic field. Recent 
MHD simulations suggest R e^ ~ 100 (|Sano fc Stone! 
I2002t lllgner fc Nelson! 120061 ) which corresponds roughly 
to an electron fraction, \ ~ 10~ 12 , although the exact 
value is dependent on the disk model adopted. Using 
our chemical model results, we calculate the value of the 
magnetic Reynolds number, ReM, everywhere and iden- 
ti fy possible dead zones where R eM ^5 100. 

IChiang fc Murrav-Clavl ((20071) argue that a second cri- 
terion, taking into account the influence of ambipolar 
diffusion in suppressing the MRI, must be applied for 
protoplanetary disks. They define a dimensionless pa- 
rameter, Am, which describes the degree to which neu- 
tral H2 molecules (which make up the bulk of the gas), 
are coupled to the accreting plasma. In order for neu- 
tral gas to be unstable, a H2 molecule must collide with 
enough ions within the e-folding time of the instability, 
l/fi, where O is the Keplerian velocity of the gas. 

Am = > Am cM (9) 

Here, x e is the disk ionization fraction, n is the number 
density of the neutral gas and /? « 2.0 x 10~ 9 cm 3 s" 1 is 
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the rate coefficie nt for ion-neutral collisi ons. Early MHD 
simulations by (|Hawlev fe Stond I1998D suggested that 
sufficient turbulence and angular momentum transport 
is achieved when the 'ambipolar diffusion parameter', 
Am, exceeds a critical value , Am cr it * 100. More recent 
simulations bv iBai fc Stond (|2011l ) suggest that in very 
weakly ionized media where the 'strong coupling' limit 
holds, such as protoplanetary disks, the critical value for 
Am is f=a 1. The 'strong coupling' limit is defined as when 
the ion inertia is negligible and the recombina tion time is 
much shorter than the orbital time (Bai 2 0111 ) . We adopt 
this latter criterion in our determination of the suscepti- 
bility of the disk to suppression of the MM by ambipolar 
diffusion. The results are reported in Section [3] 

Note, the extent and location of a dead zone is also 
dependent on the treatment of gas- grain interactions 
and the size distribu tio n of grains j San o et all l2000l : 
lllgner fe Nelsonl 120061 ). iSano et all ((20000 calculated 
that, for a fixed gas-to-dust mass ratio, the dead zone 
shrinks as the grain size increases, assuming all grains 
have the same radius. As the grains increase in size 
(likely due to coagulation) the ion density increases since 
the total surface area of grains decreases. This subse- 
quently decreases the recombination rate of gas-phase 
cations on grain surfaces. They also determine that in 
the midplane, grains are the dominant charged species 
(with charge ±e). 

It has also been postulated that gravitational grain 
settling (or sedimentation) towards the disk mid- 
plane influences the protoplanetary disk physical and 
chemical structure and thus, ionization fraction, the 
persistence of the MRI and the subs equent loca- 
tion and extent of dead zo nes (see e.g.. iChianget all 
200lUDullemond fc Dominikll2004l: lD'Alessio et al 11200a 
Nomura et al.l l2007t iFogel et all 120111 : iVasyunin et all 
20111 ). Since the dust is the dominant source of opac- 
ity in disks, grain settling allows deeper penetration of 
stellar and interstellar UV radiation potentially ionizing 
a larger prop ortion of the gaseous co mponent of the disk. 
In particular. iD'Alessio et all(|2006l ) find that an absence 
of small grains in the upper disk layer, due to sedimenta- 
tion, enhances the ionization in the disk surface and also 
decreases the temperature of the gas in the disk mid- 
plane since there is a decrease in the amount of radiation 
processed by grains and directed towards the midplane. 
In reality, disk ionization, turbulence by MRI and dust 
settling are coupled and t hus, ideally should be solved 
self-consistently (see e.g., iFromang fc Papaloizoul [20061 : 
ICieslal[2007l: iTurner et al.ll2010D . 

In this work, for our gas-grain interactions, we assume 
a constant grain radius of 0.1 /im and a fixed dust-grain 
fractional abundance of 2.2 x 10 -12 . We assume that all 
grains are negatively charged and allow the recombina- 
tion of cations on grain surfaces. This is valid assump- 
tion since negatively-charged grains dominate in regions 
where the number density, n, is < 10 12 cm and this 
holds throughout most of our disk model. We intend 
to investigate the effects of adding neutral grains and a 
variable dust-grain size distribution caused by coagula- 
tion and settling in future models. 

We have also neglected the thermal ionization of alkali 
metals, such as Na + and K + , which becomes a significant 
source of ionizati on when the gas temp e rature is greater 
than w 10 3 K (|Fromang et all l2002t lllgner fc Nelson] 



2006}). In 

our disk model, we find in the midplane the gas 
temperature reaches values higher than this only within 
a radius of a few tenths of an AU, nevertheless, it should 
looked at in future models since, due to this source of 
ionization, the gas may be magnetorotationally unstable 
close to the star. 

3. RESULTS 

We calculate the chemical abundances in the disk as 
a function of radius, height and time. The results dis- 
played in this section are those extracted at a time of 
10 6 years, the typical age of visible T Tauri stars with 
accompanying protoplanetary disks. Throughout this 
section, fractional abundance refers to the abundance of 
each species with respect to total particle number den- 
sity. As in Paper I, we ran several different models with 
differing chemical ingredients in order to determine the 
influence of each chemical process and these are listed in 
Table 01 Here, model UV-old, is our 'fiducial' model in 
which we use the same method as in Paper I to calculate 
the photochemical ra tes and we use the X -ray ionization 
rates as calculated in lNomura et all (|2007f ) . We compare 
the results from model UV-old with those from mod- 
els which include a recalculation of the photochemical 
rates only (model UV-new), the X-ray ionization rates 
only (model XR+UV-old) and both processes (model 
XR+ UV-new). All models include freeze out, thermal 
desorption, cosmic-ray-induced desorption, photodesorp- 
tion and grain-surface chemistry. 

3.1. Column Densities 

In Figure [3J we present the radial column density 
(cm -2 ) of molecules detected or searched for in pro- 
toplanetary disks at both (sub)mm and infrared wave- 
lengths. The solid lines and dashed lines are the gas- 
phase and grain-surface column densities, respectively. 

Molecules whose column densities are relatively unaf- 
fected by the method employed to calculate the photore- 
action rates and the X-ray ionization rate include CO, 
CN, CS, C 2 H, C 2 H 2 , H 2 CO, HCN, OH, CH 4 and SO. 
Over the radial extent of the disk, the column densities 
of the listed species vary, at most, by a factor of a few 
between chemical models. The gas-phase column den- 
sities of H2O and CH3OH are affected only beyond a 
radius of « 1 AU with that of CO2 altered beyond a ra- 
dius of « 10 AU. In all three cases we see a rise in the 
column density of each species when the photorates are 
recalculated. For the molecular ions, HCO + and N2H+, 
we see a different behaviour with the column density of 
the former species affected significantly within a radius 
of approximately 1 AU only. For the latter molecule, the 
column density is affected throughout the radial extent 
of the disk. For all molecules which possess an appre- 
ciable grain-surface column density, the values are rel- 
atively unaffected by the method used to compute the 
photochemical rates and the X-ray ionization rate. 

Overall, the recalculation of the photoreaction rates 
has a much bigger effect on the column densities than 
that of the X-ray ionization rate with models UV-old 
and XR+UV-old, on the whole, producing similar val- 
ues and behaviour. This is also true for models UV- 
new and XR+UV-new. Exceptions to this include N2H + 
(throughout the disk) and OH (in the outer disk beyond a 
radius of « 50 AU) . A further general observation is that 
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TABLE 4 

Chemical Models 



UV-old UV-new XR+UV-old XR+UV-new 



Chemical Process 

Thermal desorption / 
Cosmic-ray-induced desorption / 
Photodesorption / 
Grain-surface chemistry / 
Photochemistry 
X-ray ionization 

the recalculation of the photorates affects each molecule 
in a different manner. The spectrum-dependent pho- 
torates do not only directly affect the abundance and 
distribution of atoms and molecules which can undergo 
photoionization and photodissociation, they also indi- 
rectly affect the subsequent gas-phase chemistry, leading 
to enhancements/depletions of molecules which are not 
directly formed or destroyed via a photochemical route, 
e.g., HCO+. 

We see an interesting structure in the column densities 
of the sulphur-bearing species, CS and SO. Both show a 
peak between w 2 AU and ~ 10 AU. In the outer disk, 
beyond ~ 10 AU, sulphur exists primarily in the disk 
midplane as H2S ice on the grain mantle. Within 10 AU, 
the disk midplane is warm enough for H2S to evaporate 
from the grain mantle replenishing the gas with sulphur 
which forms SO and CS. Within « 2 AU, S0 2 takes 
over from SO as the dominant gas-phase sulphur-bearing 
species. CS increases in abundance again within ~ 1 AU 
in the very warm, dense midplane. 

The radial column densities, N, of a selection of 
molecules are listed in Table [3] at radii of 0.1, 1, 10 and 
100 AU. Molecules whose column densities are affected 
by more than one order of magnitude are highlighted 
in bold text. The general trend we see (with several 
exceptions discussed below) is a reduction in molecular 
column densities in the inner disk (within 1 AU) and 
an increase in the outer disk (beyond 1 AU) in models 
UV-new and XR+UV-new relative to model UV-old. At 
0.1 AU, V(HCO + ) is reduced by two orders of magni- 
tude with jV(OH) reduced by one order of magnitude, 
when the photochemistry is recalculated. The column 
density of N2H + is affected by both X-ray ionization and 
photochemistry with the former increasing the value rel- 
ative to model UV-old and the latter reducing it. At 
1 AU, only V(N2H + ) is significantly affected where it is 
enhanced by almost two orders of magnitude by the recal- 
culation of the photochemistry and reduced by a factor 
of a few by the recalculation of the X-ray ionization rate. 
At 10 AU, the molecules affected are H 2 0, C0 2 , CH 3 OH 
and N2H + . The column density of H2O and CO2 are 
both enhanced by over an order of magnitude in models 
UV-new and XR+UV-new. As at 1 AU, X-ray ionization 
reduces V(N2H + ), whereas photochemistry enhances it. 
It is a similar story for methanol, however, X-ray ioniza- 
tion has a much smaller effect than photochemistry. At 
the final radius we consider here, 100 AU, again, N2H+ 
and CH3OH are affected in a similar manner as at 10 AU. 
We also find the column density of OH increased sig- 
nificantly in model UV-new with that for CO2 also en- 
hanced in models UV-new and XR+UV-new. The only 
molecules for which X-ray ionization significantly affects 
the column density are N2H + and CH3OH. 



/ / / 

/ / / 

/ / / 

/ / / 

/ / 

/ / 



The molecular column densities, although useful for 
tracing the general radial structure and identifying those 
molecules significantly affected by the inclusion or omis- 
sion of each process, provide little information on the 
spatial distribution and abundance. This ultimately af- 
fects the strength of the line emission from the disk since 
this is influenced by the physical conditions in the region 
where each molecule is most abundant. In the following 
Section, we look more closely at the effects on the two- 
dimensional molecular structure of the disk due to the 
recalculation of the photoreaction and X-ray ionization 
rates. 

3.2. Photochemistry 

In Figure [4] we display the fractional abundance of 
HCO+, OH, H 2 0, C0 2 and N 2 H+ as a function of disk 
radius and height (scaled by the radius) for models UV- 
old (left column) and UV-new (right column). These 
molecules are those we have identified as being most af- 
fected by the recalculation of the photorates. 

In the plots for HCO + , it is clearly seen that within a 
radius of w 1 AU, the fractional abundance of HCO + in 
the 'molecular layer' in this region, located at Z/R = 0.1, 
is around two orders of magnitude lower in model UV- 
new than in model UV-old. Also clearly visible is the 
reason for the larger column density calculated beyond 
this radius, the depth of the layer of HCO + in model 
UV-new is much larger than that in model UV-old al- 
though the maximum fractional abundance attained in 
both models is similar (x(HCO + ) ~ 10 -6 ). The abun- 
dance of HCO + is controlled by ion-molecule chemistry 
and thus depends on the abundance of ionic and neutral 
precursors. An example of an ion-molecule gas-phase re- 
action which leads to the production of HCO + is 

CO + H+— ^HCO+ +H 2 . 

In model UV-new, where we see an increase in HCO + , 
we see a corresponding increase in CO and H 3 + . CO is 
directly influenced by the photochemistry since it can be 
photodissociated to produce C and O. H3" 1 ", on the other 
hand, is primarily formed via the reaction of H2 with 
H2 + and destroyed via electron recombination. The en- 
hancement in HCO + in model UV-new also corresponds 
to where we see a slight decrease in electron abundance 
(see Figure [7]). Relating this back to the photochemistry, 
this indicates an overestimation in this region in both the 
photodissociation of CO and photoionization (in general) 
in our fiducial model (model UV-old). 

The depth of the molecular layer where OH reaches its 
maximum fractional abundance (i(OH) ~ 10~ 4 ) within 
1 AU is smaller in model UV-old than in model UV- 
new, whereas, beyond this radius, the depth is larger. 
This accounts for the smaller column density in model 
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Fig. 3. — Column density (cm 2 ) as a function of radius, R, for a range of molecules detected or searched for in protoplanetary disks. 
The dashed lines represent grain-surface (ice) column densities. 
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UV-new at 1 AU and the larger at 10 and 100 AU. The 
distribution of OH is also somewhat different in model 
UV-new reaching an order of magnitude higher fractional 
abundance in the disk surface throughout the radial ex- 
tent. The abundance of OH is directly controlled by the 
photochemistry since OH is one of the products of the 
photodissociation of H2O and can itself be photodissoci- 
ated to form O and H. We see in the corresponding plots 
for water in Figure [4] that its fractional abundance is also 
slightly higher in the disk surface in model UV-new rel- 
ative to model UV-old. We also find that the abundance 
of free oxygen atoms decreases in this region, indicating 
in model UV-new, more atomic oxygen is locked up in 
oxygen-bearing molecules than in model UV-old. 

Looking at the distribution of water, we see a large en- 
hancement in the fractional abundance in the molecular 
layer beyond a radius of w 1 AU, going from a value of 
~ 10~ 6 in model UV-old to ~ 10~ 4 in model UV-new. 
Not only is the maximum fractional abundance reached 
much higher, the extent over which water exists with a 
value > 10 -9 in model UV-new is also much larger. We 
conclude that the abundance of gas-phase water is very 
sensitive to the method employed to calculate the pho- 
torates. 

Both OH and H2O can be formed via neutral- neutral 
gas-phase r eactions in warm reg ions of the disk where 
T > 200 K UGlasseold et al.ll2009h . 

H 2 + — >OH + H 
H2 + OH — >H 2 + H 

Since we see an increase in the fractional abundance of 
H2 in model UV-new relative to UV-old over the region 
where both OH and H2O are increased, this gas-phase 
production route for both species is more important in 
model UV-new than in model UV-old. Hence, the en- 
hancement seen in the abundances of OH and H2O in 
model UV-new is a combination of increased gas-phase 
production and decreased photodestruction. In Figure El 
we display the photodissociation rates of OH and H2O as 
a function of disk height (scaled by the radius) at radii 
of 1 AU, 10 AU and 100 AU. A decrease of around an 
order of magnitude in the photodissociation rates of both 
species is clearly seen in the upper disk layers account- 
ing for the increase in abundance of both species in this 
region. The gross overestimation of the the photodisso- 
ciation rates in model UV-old is due to the inclusion of 
Lyman-a photons in the calculation of the wavelength- 
integrated UV photon flux (see Section 12. 2. lj) . Also, note 
that the dissociation rates vary differently as a function 
of height in model UV-new versus UV-old since in the 
former model, the wavelength dependence of both the 
photo cross sections and UV field are included. 

Gas-phase CO2 is affected mainly in the outer disk be- 
yond a radius of ps 10 AU. We see that CO2 in this region 
exists in a layer lower than that of water (due to its lower 
binding energy to dust grains) and is enhanced in the 
outer disk from a fractional abundance ~ 10~ 6 in model 
UV-old to ~ 10 -4 in model UV-new. In model UV- 
new, in this region, it possesses a comparable fractional 
abundance to the other main oxygen-bearing species, CO 
and O2. As in the case for water, the extent over which 
2KCO2) has a value > 10~ 9 is much larger in model UV- 
new than in model UV-old. We see an enhancement in 



CO2 when the fractional abundance of ionic carbon be- 
gins to increase relative to atomiccarbon. Note, CO2 is 
destroyed via reaction with C + to form CO + and CO. In 
model UV-new, the boundary where C + /C ~ 1 is much 
higher in the disk than in model UV-old (at a height of 
w 5 AU and w 2 AU, respectively, at a radius of 10 AU). 
Thus in model UV-new, gas-phase C0 2 can evaporate 
from the grain surface and remain intact in this layer 
due to the lack of sufficient gas-phase destruction mech- 
anisms. We display the photodissociation rate of C0 2 as 
a function of disk height (scaled by the radius) at radii 
of 1 AU, 10 AU and 100 AU in Figure El As in the case 
for OH and H2O, we see an decrease of around an order 
of magnitude in model UV-new compared with UV-old. 
This accounts for the increase in the abundance of C0 2 in 
the upper disk layers in model UV-new. Note, in model 
UV-old, the dissociation rates for CO2 are consistently 
larger than that of H2O and OH, reflecting the the rates 
calculated for the unshielded interstellar medium (14, 5.9 
and 3.5 x 10~ 10 s _1 , respectively). However, in model 
UV-new, the H2O photodissociation has the largest rate 
at each radius, with OH also having a larger rate than 
CO2 at 100 AU. This demonstrates how the relative pho- 
torates are sensitive to the shape of the radiation field 
due to the specific variation in the photo cross section of 
each species. 

When the photochemistry is recalculated, we see a dra- 
matic increase in the region over which N2H + reaches its 
maximum fractional abundance, now residing in a molec- 
ular layer which permeates the entire disk. The most 
dramatic increase is seen in the inner disk (< 100 AU) 
where ir(N2H + ) jumps from a maximum value of ~ 10~ 12 
in model UV-old to ~ 10~ 9 in model UV-new. Similar 
to HCO + , the abundance of N2H" 1 " is controlled by ion- 
molecule chemistry, it is formed via the reaction of N2 
with H 3 + and is destroyed effectively by reaction with 
CO and electron recombination. Looking at Figure [7] 
which displays the electron fractional abundance for all 
four models as a function of disk radius and height, we 
see in the layer where N2H + is significantly enhanced, 
there is a respective decrease in the electron fraction in 
model UV-new. N 2 H + and HCO + show different be- 
haviours since the abundance of N2H + is more sensitive 
to the abundance of H 3 + than HCO+. HCO+ has various 
alternative routes to formation other than the reaction 
of CO and H 3 + , e.g., C + + H 2 0, whereas N 2 H+ forms 
mainly via the reaction of N2 with H 3 + . We discuss 
N2H + chemistry further in Section [3~3l 

In Figure El we present the fractional abundances of 
several atoms and ions at a radius of 1 (top panel) and 
10 AU (bottom panel). We display the results for C + 
also since the ionization potential for carbon is much 
lower than that for oxygen and nitrogen. In model UV- 
new, there is a significant depletion of free O atoms and 
C + ions in the region where we see an increase in the 
fractional abundance of OH, H2O and the main carbon 
species, C. Note that the atomic carbon abundance in 
the disk surface is larger in model UV-new than in model 
UV-old, however, here free carbon exists mainly in ionic 
form as C + . In this region, we also see an enhance- 
ment in the fractional abundances of the molecular ions, 
HCO + and N2H" 1 ". Hence, in model UV-new, in the 
molecular layer where molecules reach their maximum 
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fractional abundance, more oxygen atoms are contained 
within molecules than in model UV-old. Also, in this 
region, less carbon exists as C + so the ratio of C + /C is 
generally lower in model UV-new than in model UV-old. 
This can be attributed to (1) the photorates in model 
UV-old being overestimated for particular molecules and 
atoms and (2) an adjustment of the gas-phase chemistry 
due to the enhancement / depletion of species directly af- 
fected by the radiation field. 

The layer in which there is an enhancement in the 
abundance of molecules also coincides with the boundary 
where we cross from hydrogen in mainly molecular form 
(as H2) to atomic form. In model UV-new, this transi- 
tion is smoother and results in a slightly lower fractional 
abundance of H atoms over the extent of the molecular 
layer. At high temperatures, the destruction of molecules 
by free hydrogen atoms can become an important pro- 
cess. The region where we find the large increase in OH 
and H2O, coincides with a gas temperature of > 1000 K 
so that the reduction of free H atoms (which act to de- 
stroy molecules) in model UV-new is also influencing the 
gas-phase chemistry and subsequent abundances. 

3.3. X-ray Ionization 

We find the only molecule significantly affected by X- 
ray ionization is N2H" 1 ". We see a decrease in both the 
maximum fractional abundance attained by N2H + and a 
reduction in the spatial extent over which it exists with 
an appreciable abundance. Note, in model UV-old, the 
distribution of N2H 4 " is also different to that of any of 
the other molecules discussed thus far, existing mainly 
in the outer disk beyond a radius of 100 AU and ex- 
tending into the disk midplane (see Figure [4]) . Since the 
abundance of N2H+ is controlled by ion-molecule chem- 
istry, the relative abundances of precursor ions and neu- 
tral molecules will influence the amount of N2H 4 " which 
can exist. N2H + is formed in cold, dense regions via 
reaction of H3 + with N2, the latter of which is formed 
via radical-radical reactions e.g., N + NH — > N2 + H. 
The abundance of H3 + is controlled by the ionization of 
H2 to form H2 + . In reactions where molecular ions can 
form via proton donation, the ionization rate can have 
a large impact on the resulting abundances reached, as 
seen here for the case of N2H+. Also, N2H 4 " is effectively 
destroyed by electron recombination so that the ioniza- 
tion fraction in the disk plays an important role (see Sec- 
tion [33J. However, we do not see this effect in the abun- 
dance and distribution of HCO + which also depends, al- 
beit to a lesser extent (see later), on the abundance of 
H3 + and which appears independent of the treatment of 
X-ray ionization. HCO + resides in a different layer in 
the disk where the abundance and distribution of H3 + 
and electrons do not vary significantly between models 
UV-old and XR+UV-old. We do see corresponding dif- 
ferences in the abundance and distribution of H 3 + , in 
particular, between models UV-old and XR+UV-old in 
the region where N2H + is most abundant. HCO + also 
has a plethora of different routes to formation, other than 
the protonation of CO by H3 + . 

In our disk model, X-ray ionization rate dominates in 
the outer disk upper layers over cosmic-ray ionization. 
Our recalculation of the rates in this region leads to a 
decrease in N2H + indicating an overestimation of the X- 
ray ionization rate in our original model. 



3.4. Disk Ionization Fraction 

In Figure[7]we display the fractional abundance of elec- 
trons as a function of disk radius and height for all mod- 
els considered in this work and in Figure [51 we present 
the fractional abundances of several abundant cations as 
a function of disk height at radii of 15 AU (top) and 
150 AU (bottom). 

At first glance, there appears little difference in the dis- 
tribution of electrons between the different models. The 
ionization fraction in the disk ranges from close to 1, in 
the disk surface nearest the star, to less than 10 -12 in the 
disk midplane between a radius of 1 AU and 10 AU. Com- 
paring the plots for models UV-old and XR+UV-old, we 
see only minor differences; in the boundary layer where 
the ionization gradient is highest and coinciding with 
the aforementioned molecular layer, and in the midplane 
layer where the electron fraction reaches values < 10 -11 
which is slightly larger in depth in model XR+UV-old. 

Comparing models UV-old and UV-new we see there is 
more of an effect in the boundary layer described above 
which extends into the outer disk, in fact, we see that 
at the outermost radius, ~ 300 AU, the ionization frac- 
tion in the disk surface layer (Z/R rs 0.6 - 0.7) in model 
UV-new is orders of magnitude lower than that in model 
UV-old. This is coincident with where we see an increase 
in the abundance of N2H" 1 ". We also note the slight in- 
crease in the electron abundance in the disk surface in 
model UV-new compared with model UV-old. This is 
unsurprising, as we would expect the stellar radiation 
field to impact on the chemistry in the disk surface, in 
particular. Comparing the graphs for models UV-new 
and XR+UV-new, we see the effects of the X-ray ioniza- 
tion in the disk surface, which negates the decrease in 
the electron abundance to some extent. 

Looking at Figure [FJ in the disk surface at 15 AU, 
the most abundant cations are H + and He + whereas at 
150 AU, the most abundant is C + . In the disk midplane, 
at all radii > 1 AU, HCO + dominates the cation frac- 
tion. Within ~ 1 AU, metallic ions such as Fe + and Na + 
also contribute to the cationic abundance in the midplane 
forming efficiently via charge exchange (at the expense 
of molecular cations such as HCO + ). In the molecular 
layer, the picture is more complex with several molecu- 
lar and atomic ions contributing to the cation fraction 
including H 3 +, H 3 + , HCO + and C + . In this region, 
at 150 AU, metallic ions with an appreciable fractional 
abundance (~ 10 -9 to 10~ 8 ) such as Fe + , Mg + and Na + 
also contribute to the cation fraction. N2H+ only reaches 
significant abundances in the the outer disk (> 100 AU) 
caused, to a degree, by the depletion of gas-phase CO 
due to freeze out onto grain surfaces in this region. 

Following the theory outlined in Section 12. 3\ we cal- 
culated the magnetic Reynolds number, Rcm, and the 
ambipolar diffusion parameter, Am, as a function of disk 
radius and height and our results are displayed in Fig- 
ure 02 We plot the results for model XR+UV-new only, 
since we find both the photochemistry and X-ray ioniza- 
tion have little impact on the values of Reyi and Am in 
the inner disk mid plane. 

Simulations by iSano fc Stond (|2002f ) suggest that a 
suitable value for Re^ is ~ 100. At values lower than 
this, accretion in the disk is likely inhibited due to sup- 
pression of the MRI by Ohmic dissipation. Looking at 
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Fig. 5. — The photodissociation rates of OH, H2O and CO2 as 
a function of disk height (scaled by the radius i.e., Z/R) at radii 
of 1 AU, 10 AU and 100 AU for models UV-old (solid lines) and 
UV-new (dashed lines). 
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Fig. 6. — The fractional abundance of C + , C, O and N as a 
function of height, Z, at radii of 1 AU (top) and 10 AU (bottom) 
for models UV-old (solid lines) and UV-new (dashed lines). 



the left panel of Figure O and adopting this criterion, we 
see that there is a probable dead zone where the MPJ 
is susceptible to Ohmic dissipation in the midplane ex- 
tending from the innermost radius, 0.1 AU, to w 20 AU. 
The depth of this dead zone varies from Z/R = 0.1 at 
0.1 AU to Z/R w 0.15 at 10 AU. We find little differ- 
ence in the extent of the dead zone between the various 
chemical models. 

Similarly, the simulations of iBai fc Stone! (|2011[ ) sug- 
gest that ambipolar diffusion can inhibit mass and an- 
gular momentum transport in regions in the disk where 
the ambipolar diffusion parameter, Am, exceeds a crit- 
ical value, Am crlt ~ 1. In the right-hand panel of Fig- 
ure [9l the white contours represent the boundary where 
Am = 1 whereas the black contours represent Am = 100. 
We see ambipolar diffusion has a bigger effect creating 
a much larger dead zone in the disk midplane, ranging 
from just beyond 2 AU ou t to around 200 AU. I f we 
adopt the earlier criterion of lHawlev fe Stonel (|1998|) , i.e. 
Am clit = 100 (black contours), we find that the MRI is 
s uppressed throughout most of the disk mass. 

iPerez-Becker fe Chiand (|2011[ ) argue that, in order for 
efficient accretion in protoplanetary disks, both criterion 
for ReM and Am must be met. Our calculations sug- 
gest, in our disk model, most of the midplane is inactive 
and so accretion can occur via t he surfac e layer s only, 
a similar result to that found by iGammid (fl99 6). As a 
measure of the effect on the mass accretion rate of the 
disk due to the presence of dead zones, in Figure [10] wc 
show the ratio of the 'active' column density to the to- 
tal column density, Sact/^tot, as a function of disk 
radius, R, since this ratio estimates the efficiency of ac- 
cretion. We plot Sact for four criteria: ReM > 100 (solid 
red line), Am > 100 (solid green line), Am > 1 (dashed 
green line), and ReM > 100 and Am > 1 (dotted blue 
line). For the latter criterion, we see that accretion effi- 
ciency ranges from < 1 % in the inner disk (< 10 AU) 
to around 70 % in the outer disk (> 200 AU). Although 
most of the outer disk is unaffected by Ohmic dissipation 
(solid red line), it is susceptible to ambipolar diffusion 
(dashed green line). Applying the more stringent condi- 
tion of Am > 100 (solid green line), we find that most 
of the disk is constrained to an accretion rate < 1% that 
which dictates the disk structure. The structure seen 
in the active column density is due to the corresponding 
structure seen in the spatial distribution of ReM and Am 
(displayed in Figure HJ which are related to the fractional 
and absolute electron density, respectively. 

Compared with Ohmic dissipation, ambipolar diffusion 
suppres ses the MRI in lower d ensity regions (see e.g., Fig- 
ure 1 in lKunz fe Balbusll2004f) . Our results are consistent 
with this since we see the dead zone extending into the 
lower density regions in our disk i.e., to larger radii in 
the disk midplane and towards the disk surface, when 
the ambipolar diffusion criterion is included. 

The suppression of accretion in specific areas leading to 
a build up of material flowing through the disk midplane, 
will change the physical nature of the disk. This, in turn, 
influences the chemical composition and thus the pen- 
etration of UV photons, X-ray photons and cosmic-ray 
particles. The changing physical and chemical conditions 
may act to increase ReM and Am to the extent that they 
exceed their critical values, hence, the presence of dead 
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Fig. 7. — The fractional abundance of electrons as a function of disk radius, R, and height (scaled by the radius i.e., Z/R) for models 
UV-old (top left), UV-new (top right), XR+UV-old (bottom left) and XR+UV-new (bottom right). 



zones may be time dependent. Young stars and protostel- 
lar objects often undergo transient periods of enhance- 
ments in luminosity often ascribed to spells where the ac- 
cretion disk possesses an increased mass flow rate going 
from < 10 ~ 7 M Q yr" 1 to ~ 10 ~ 4 M yr" 1 for periods 
of ~ 100 yr. An episodic build up and release of material 
in the disk due to the development and subsequent dis- 
sipation of an accompanying dead zone may be respon- 
sible (see e.g. . iHerbid 119771; lHartmann fc Kenvonl 119961: 



Calvet et aljl2000HArmitage et al.ll2001l : iVorobvov et al l 
200G). 

We note here that the steady disk model is a first step 
in our attempt to approach reality. A more sophisticated 
model, in which we treat the charge balance more real- 
istically (including detailed balance of grain charge and 
PAHs) and in which we calculate the dust evolution in 
conjunction with the chemical evolution, may go some 
way to addressing the apparent discrepancy between the 
observed and calculated accretion rates in the inner disk. 
We also note here, the disk mass and surface density pro- 
file of our model is consistent with observed thermal dust 
emission from the outer regions of protoplanetary disks 
and observed accretion rates are mainly derived from op- 
tical and UV spectral li nes originating fro m close to the 
stellar surface (see e.g.. lCalvet et al.1l2000l ). 

The suppression of turbulence (and thus angular mo- 
mentum transport) has various further consequences on 
the disk physical structure and resulting chemistry. A 
significant source of heating in the disk midplane close 
to the star is viscous dissipation due to accretion. This 
can elevate the temperature in this region to ~ 1000 K 



which has implications on the disk chemical structure 
and resulting line emission. Since this is typically the 
densest region of the disk, a lot of material is contained 
within a few AU of the star and thus can have a signifi- 
cant contribution to the d isk-integrated line profiles (see 
e.g.. iCarr fc Najital 120081 ) . If accretion is suppressed in 
this region, this heating cannot occur and so will affect 
the chemistry and resulting line emission from the inner 
disk. It is possible that there are alternative sources of 
ionization in the very hot, dense inner regions of disks, 
such as thermal ionization, which allow the gas to be- 
come ma gnetorotationally unstable ver y close to the star 
(see e.g.. lUmebava shi fc Na kano 1988|). 

We have investigated t he effect of turbulent m ixing on 
disk chemical structure (|Heinzeller et al.ll2011l) and find 
that in the 'planet-forming' region, within approximately 
10 AU and where the bulk of the infrared line emission 
originates, turbulent mixing acts to mix material in the 
vertical direction, dredging material from the disk mid- 
plane towards the surface. In fact, we find that our mod- 
els with turbulent mixing included agree best with cur- 
rent infrared data from disks. In dead zones, the disk flow 
is essentially laminar and so turbulent mixing within the 
dead zone itself may be suppressed. This could act to 
confine a proportion of infra-red active molecules to the 
cold, dense midplane. 

Finally, the only negative charge carriers we consider 
here are electrons and simple negative ions such as H" 
and CN~. Dust grains and polycyclic aromatic hydro- 
carbons (PAHs) are also able to capture electrons to form 
negatively charged grains and PAHs. As discussed pre- 
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Prairie k Bertoldi I 119961: iLee et all 119961: iVisser et all 
2009). In protoplanetary disks, detailed two-dimensional 
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Fig. 8. — The fractional abundance of several important cations 
as a function of height, Z, at radii of 15 AU (top) and 150 AU 
(bottom) for model UV-new. 



viously in Section |2~5I dust grains may be the dominant 
charge carriers in the disk midplane. Effective accretion 
depends on the degree of coupling of charged particles to 
the magnetic field. If PAHs and dust grains are the dom- 
inant charge carriers, the relatively large, heavy PAHs 
and dust grains will increase drag, further opposing ac- 
cretion. The addition of charged grains is a modification 
we intend to investigate in future models. 

3.5. Self and Mutual Shielding of Hi and CO 

The photodissociation of H2 and CO is dominated by 
line absorption rath er than continuum absorption (see 
e.g., ILee et afl 119961 ). hence, in any atmosphere where 
there is a significant column density of foreground H2 or 
CO, there will be a degree of self shielding against pho- 
todissociation dependent on the amount of intervening 
material. Also, there is an overlap in wavelength ranges 
over which both molecules absorb leading to a degree of 
mutual shielding due to the removal of UV flux by the 
photodissociation of foreground H2. 

To date, significant investigation into the self- and 
mutual shielding of H2 and CO in photon-dominated 
regions (PDRs) has been conducted leading to a va- 
riety of approximate methods for the computation of 
self-shielding factors in plane-parallel media (see e.g., 



radiative transfer is required since the radiation field 
has two sources, the central star and the interstellar 
medium. Also, the addition of self-shielding requires that 
the abundances of H2 and CO are known everywhere in 
the disk at each time step in the calculation simultane- 
ously and this is computationally impossible in a spa- 
tially high-resolution model such as this. A compromise 
often employed in disks is to assume the plane-parallel 
approximation in the vertical direction and calculate the 
chemical structure from the disk surface downwards. In 
this way, the column density of H2 and CO from the disk 
surface to the point of interest can be estimated by ex- 
tracting out the abundances in the upper layers at a par- 
ticular time (e.g., 10 6 years). The self-shielding factors 
calculated for plane-parallel media are then employed in 
the calculation of the photodissociation rates. However, 
in the upper disk, where photodestruction is most sig- 
nificant, the direct stellar radiation field is the dominant 
component so that the column density of material from 
the star to the point of interest becomes the important 
parameter opposed to t he column density from the disk 
surface (see Figure 4 in INomura k Mmarl (l20ul) 1. Also, 
many self-shielding factors are computed assuming irra- 
diation by the interstellar field only. As we have already 
discussed, the radiation field in protoplanetary disks will 
have a very different shape and strength to that found in 
irradiated interstellar clouds. The computation of such 
factors also depends on the assumed gas temperature and 
density since these parameters affect the line shapes and 
thus absorption cross sections of H2 and CO. In sum- 
mary, the treatment of self- and mutual shielding in pro- 
toplanetary disks is a complex problem and the use of ex- 
isting approximations to model the self-shielding in disks 
is problematic as the application of shielding factors com- 
puted for irradiated interstellar clouds is debatable. 

In the generation of our disk phy si cal m odel, 
which originates from INomura k Millar] (|2005), the 
self-shielding of H 2 is included in the calculation 
(|Federman et al. ill 9 791 ). Here, althou gh we use the cal- 
culate d H and H2 abundances from INomura k Millar 
(2005) as initial conditions in our chemical model, the 
abundances of both species are adjusted to an extent 
when the full, and vastly more complex, time-dependent 
chemistry is computed. 

In order to quantify the effect that self- and mutual 
shielding of CO may have in o u r mod el, we adopt the 
method described in ILee et all (| 19961 ) to compute the 
photodissociation rate of CO as a function of disk height 
at a radii of 1 AU, 10 AU and 100 AU assuming the 
usual plane-parallel approximation. The results of this 
calculation are shown in Figure [TTJ The shielding factors 
only become significant deep in the atmosphere where 
the radiation field decreases significantly in strength. We 
conclude, the addition of self-shielding of CO using ap- 
proximation methods has a minor effect on the overall 
disk structure when compared with the effects of the re- 
calculation of the photorates taking into consideration 
the two-dimensional radiation field in the disk. We ex- 
pect the inclusion of self- and mutual shielding in the 
radial and vertical directions will lead to an increase in 
the abundances of H2 and CO in the outer disk, in partic- 
ular, since the radial column densities here will be larger 
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Fig. 9. — The magnetic Reynolds number, Reyi (left), and ambipolar diffusion parameter, Am (right), as a function of disk radius, R, 
and height (scaled by the radius i.e., Z/R) for model XR+UV-new. The black contours on both plots represent the boundary where Re^ 
and Am = 100, whereas, the white contours on the right-hand plot represent the boundary where Am =1. 
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Fig. 10. — The ratio of the active column density to the total 
column density, Eact/£tot, as a function of radius, R, for four 
criterion: Re^ > 100 (solid red line), Am > 100 (solid green line), 
Am > 1 (dashed green line), and Reyi > 100 and Am > 1. 

than the vertical and the radial stellar radiation field is 
dominant. We intend to look more closely at the com- 
plex issue of self- and mutual shielding of H2 and CO in 
disks in future work. 

3.6. On the Effects of Lyman-a Radiation 

In Section [2. 2 .11 we discussed how we included Lyman- 
a radiation in the calculation of the wavelength- 
integrated flux yet neglected it in the determination of 
the wavelength-dependent UV spectrum throughout the 
disk. This was due to the difficulty in treating Lyman-a 
scattering in protopl anetary disks, an issue on ly recently 
addressed in work bv lBethell fc Berginl ()2011[ ). Here, we 
consider the effects of the addition of Lyman-a on the 
calculated photodissociation rates of several molecules. 

We calculate the photodissociation rates using the stel- 
lar spectrum displayed in the left-hand panel of Figure [1] 
which is the unshielded UV photon flux at the disk sur- 
face at a radius of 1 AU. Note, this spectrum is different 
to that used to calculate the photorates displayed in Fig- 
ure [5] In our full disk model, we assume that the stellar 
radiation has travelled through, and thus been attenu- 
ated by, material along the line of sight between the star 
and the disk surface. 
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Fig. 11. — CO shielding factors (solid lines) and photodisso- 
ciation rates (s^ 1 ) with (dot-dashed lines) and without (dotted 
lines) the shielding factors implemented as functions of disk height 
(scaled by the radius i.e., Z/R) at radii of 1 AU, 10 AU and 100 AU. 

We calculate the unshielded photorates with and with- 
out Lyman-a, using t he calculated Lyman-a cro ss sec- 
tions from Table 1 in van Dishoeck etaTl (|2006D . Our 
results are presented in Table [SJ In our determination 
of the photorate at the Lyman-a wavelength, we assume 
the cross section is constant across the line, determine 
the rate due to Lyman-a only, and add that to the rate 
calculated by integrating over the background UV field 
(see Equations[3]and|4|). The integrated photon flux over 
the Lyman-a line is 8.4 x 10 14 photons cm~ 2 s _1 . 

Note that the photodissociation cross section of most 
molecules at the Lyman-a wavelength is generally lower 
than at other wavelengths, hence, we find that the inclu- 
sion of Lyman-a increases the photodissociation rates by 
no more than a factor of 2.52 (as determined for HCN). 
In general, the enhancement is less than a factor of 2 
which is well within the uncertainties in both the UV 
spectrum and the photo cross sections. Even though 
the UV field is weaker at wavelengths other than the 
Lyman-a wavelength (see Figure [T]) the overall contribu- 
tion of Lyman-a photons to the photodissociation rates 
is around the same order of magnitude as the contri- 
bution by the background UV photons, i.e., we do not 
see the order of magnitude difference in the photorates 
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when the Lyman-a photon flux is included in the inte- 
grated UV field and the interstellar rates scaled by this 
flux (se e Figure |5l). This r esult is in agreement with the 
work of iFogel et al.1 (|2011h who find that the inclusion of 
Lyman-a has an effect on the chemistry only when grain 
settling is included in their model. 

3.7. Comparison with Other Models 

A full two-dimensional calculatio n of the UV radia- 
tive tr ansfer in disks was achieved bv lvan Z adclhof f et al.l 
(2003). They coupled the UV radiation field and photo- 
chemistry and applied their chemical model to the es- 
tablish ed protoplanetary disk model of iD'Alessio et all 
(1999). In this earlier work, they calculate the pho- 
torates adopting difference shapes for the stellar UV field 
and come to a similar conclusion as ourselves i.e., the 
abundances of molecules are sensitive to the shape of 
the adopted stellar UV spectrum. They also find at a 
radius of 105 AU in their model, that the photodisso- 
ciation rates of C2H and H2CO are generally underesti- 
mated when the UV stellar spectrum is approximated by 
a scaling of the interstellar radiation field. 

Early investigations into the calcul ation of the X-ray 
ionization in dis ks include the work of | Aikawa fc Herbstl 
(119991 1200 ID . In lAikawa fe Herbstl (|1999ft . thev adopt a 
power-law approximation for the X-ray a bsorption cross- 
sectio n si milar to that used in Pa per I (Malonc v et al.l 
119961 ). In lAikawa fe Herbstl (|2001l ) they perform a cal- 
culation similar to that adopted here, where the ex- 
plicit elemental composition and absorption cross-section 
are taken into account as is the direct X-ray ionization 
of elements. Although their underlying disk physical 
model is analytical in nature and different to ours in 
that the X-rays are included in purely a chemical sense 
i.e., X-ray heating of the disk is not taken into account, 
they too find little difference between results using ei- 
ther method. Mo r e rece nt mo dels based on the work o f 
lAikawa fc HCTbstl (|200l (e.g.. IWoods fe Willacvl[2009l) . 
have adopted this method of calculating the X-ray ion- 
iz ation rates. 

lAikawa fe Nom ura (2006|) include spectrum-dependent 
photochemistry looking at the effect of grain growth on 
disk chemical structure although they neglect the effects 
of X-rays. The p hysical model used i n their work is based 
on the model of (Nomura fc Millarl (|2005l) where a full 
two-dimensional calculation of the radiative transfer al- 
lowed the determination of the UV spectrum (including 
UV excess emission from the star) everywhere in the disk. 
We currently have in preparation a paper where we look 
at the effects of grain growth in a protoplanetary disk in 
which X-rays are included (Walsh et al. 2011, in prepa- 
ration) . 

The disk model of iWoitke et~aT1 (|2009ft . 'ProDiMo', 
has recently be en updated with spectrum-dependent 
photo chemistry (iKamp et al.l 120101) and X-ray ioniza- 
tion (|Aresu et all I2011D with the methods used and 
cross-sectional data adopted simila r to th at used here. 
The recent work bv lVasvunin et all (|2011[ ) now includes 
wavelength-depen dent photochemistry, updating their 
previ ous model (jSemenov et al.l 120051 : iVasvunin et all 
2008) to investigate t he effects of grain- settling on the 
disk chemistry. Also, IFogel et ail (|2011[ ) have included 
photochemistry investigating the effect of Lyman-a ra- 
diation on the disk chemical structure They find, in mod- 



els in which dust settling is included, those species with 
significant photodissociation cross sections around the 
Lyman-a line are particularly affected, generally lead- 
ing to a reduction in column density for species such as 
HCN, NH3 and CH4. This effect is significantly less in 
models where they assume the dust is well-mixed with 
the gas, an assumption we adopt in our work. 

3.8. Comparison with Observation 

The observation of molecular line emission from proto- 
planetary disks remains a challenge due to the small an- 
gular size of these objects and the limitations of existing 
(sub)mm observing facilities. Ideally, to directly com- 
pare our results with observations, full radiative transfer 
should be performed and disk-averaged line profiles and 
intensities calculated to compare with those measured by 
telescopes. 

Our model is not specific to any particular source 
and we have adopt ed the star-disk parameters of a typ- 
ical T Tauri star ( Ha rtmann fc Kenvonl [L996T ) . How- 
ever, since we adopt the X-ray spectrum and UV ex- 
cess observed from the T Tauri star, TW Hya, we com- 
pare our calculated column densities with those derived 
from molecular line observations (single dish and inter- 
ferometric) of this source. TW Hya is a well-studied, 
relatively old (~ 10 Myr), nearby (51 ± 4 pc), almost 
face-on (< 10°), iso l ated classical T Tauri star (see e.g. , 
Kastner et al.lll997t IWebb et al.|[l99l IKrist et al1[2000l: 
Mamaiekl 120051 ). Observations of continuum dust emis- 



sion from TW Hya suggest it has an cleared inner hole, 
possibly caused by a forming plane t, and a truncated 
outer disk with rad i us » 200 AU dCalvet et al.l 120021: 
iSetiawan et al.ll2008HAkeson et al.ll2011l) . 

Rotational line emission at (sub) mm wavelengths from 
several molecules has been observed in TW Hya in- 
cluding CO, CN, HCO+, HCN, as well as the isotopo- 
logue s , 13 CO, H 13 CO+, DCO+ an d DCN (IKastner etall 



19971 Ivan Zadelhoff et all 120011: Ivan Dishoeck) 
Wil ner et all 120031: lOi et al.l 120041 : iThi et al.l 



2003 



2004; 



Qi et all I200fil l2008h . Molecular line emission has been 



spatially mapped by iQi et all (|2004 120061 12008ft using 
the S ubmillimeter Array (SMA) and by I Wilner et all 
(2003) using the Aust ralia Telescope Compact Array 
(ATCA). More recently, iHogerheiide et~aTl (|2011[) report 
the detection of the ground-state emission lines of both 
spin isomers of H2 O in TW Hya observed using the Het- 
erodyne Instrument for the Far-Infrared (HIFI) on the 
Hers chel Space Observa tory. In addition to (sub)mm 
data. lNaiita et all (|2010T ) present detections of OH, CO2, 
HCO + , and possibly CH3, at mid- infrared wavelengths 
using the Infrared Spectrograph (IRS) on board the 
S pitzer Spac e Telescope. 

IThi et all (|2004h and IQi et all (|2008l ) each use a 
method to derive the column den sity of eac h mole cule 
using their observed line profiles. IThi et all (|2004f ) as- 
sume the emitting region has a typical density and tem- 
perature, constrained by line ratios in a single molecule, 
and assuming the energy levels are thermalised (i.e. they 
assume local thermodynamic equilibrium or LTE). The 
optical depth can be estimated using line ratios of the 
same transition in two isotopologues assuming the same 
excitation temperature. They list their beam-averaged 
column densities for TW Hya assuming a disk radius of 
165 AU and an excitation temperature of 25 K. We re- 
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TABLE 5 

Unshielded Photodissociation Rates at 1 AU 



Species 


CLy 1 


fc L y 


fcuv 


khy + fcuv 


(fc L y+fc UV )/fc U V 


co 2 


6.1 (-20) 


5.12 (-05) 


9.69 (-03) 


9.74 (-03) 


1.01 


H 2 


1.2 (-17) 


1.01 (-02) 


1.41 (-02) 


2.42 (-02) 


1.72 


HCN 


3.0 (-17) 


2.52 (-02) 


1.78 (-02) 


4.30 (-02) 


2.52 


OH 


1.8 (-18) 


1.51 (-03) 


6.33 (-03) 


7.84 (-03) 


1.51 


C 2 H 2 


4.0 (-17) 


3.36 (-02) 


5.16 (-02) 


8.52 (-02) 


1.65 


CH 4 


1.8 (-17) 


1.51 (-02) 


1.44 (-02) 


2.95 (-02) 


1.51 


NH 3 


1.0 (-17) 


8.40 (-03) 


2.56 (-02) 


3.40 (-02) 


1.33 



References. — (1) Ivan D ishocck ct al. (2006) 
Note. — a(b) means a x 10 h 



produ ce their calculated values here in Table [SJ iQi et al.l 
(2008) constrain their estimated column densities be- 
tween 10 AU and 100 AU by assuming a simple power-law 
distribution for the surface density of the disk and fit- 
ting the vertical extent and abundance of each molecule 
to match the observed line profiles using a x 2 method. 
Again, we list their fitted values at a radius of « 100 AU 
here in Table [6] 

Comparing our calculated column densities at a radius 
of 100 AU (listed in Table [3]) with those derived from 
observation, we see our value for CO is more than two 
orders of magnitude larger than the value from lThi et al.l 
(2004). Comparing the values for HCO+ and HCN, we 
see good agreement between the column densities from 
both authors, however, our calculated column densities 
for both species are around one order of magnitude larger 
than the observed values. Looking at CN, we see the 
same behaviour with our model predicting a value be- 
tween 2 to 5 x 10 14 cm 2 compared with the observed 
value of 6.6 x 10 13 cm 2 . 

Unfortunately, iHogerheiide et al] ()201lD . do not re- 
port their results in terms of column density, however, 
they do interpret their resul ts using the detailed disk 
model of Woi tke et al.l (|2009[ ) calculating a total water 
vapour mass of 7.3 x 10 24 g in a disk with total mass 
1.9 x 10 -2 Mq. This corresponds to a disk-averaged wa- 
ter vapour fractional abundance of ps 2 x 10~ 8 . Taking 
our modelled column densities of H2 and H2O at 100 AU 
for model XR+UV-new, we find a column-averaged frac- 
tional abundance of 7 x 1 0~ 8 which compare s well with 
the disk-averaged value. rHogcrhcij de et al.l (|2011| ) re- 
port that their observations are optically thin hence the 
line emission is tracing the entire column density of wa- 
ter throughout the vertical extent of the disk. Also, 
iNaiita et al.l (|2010l ) do not convert their line observa- 
tions to physical quantities such as column density or 
fractional abundance since they plan to disuss their re- 
sults for OH, HCO + and CO2 in more detail in a future 
publication. They do mention their apparent anoma- 
lous result for TW Hya which is their non-detection of 
H2O, C2H2 and HCN line emission, in contrast to re- 
sults for most classical T Tauri stars ob served using IRS 
(jCarr fc Najital 120111: iSalvk et al.l 120111 ) . One explana- 
tion they give is that TW Hya is unusual in that is ap- 
parently at a more advanced stage than its accretion rate 
suggests i.e., TW Hya is not an archetypical classical 
T Tauri star. 

Why is our model overestimating the column densities, 
especially those observed with ground-based facilities? 



TABLE 6 

TW Hya Observed Column Densities 



Species 



N (cm -2 ) 



Beam-averaged values from Thi et al. 



CO 
HCO+ 
CN 
HCN 



3.2 (16) 
4.4 (12) 
6.6 (13) 
9.2 (12) 



Fitted values at 100 AU from Qi et al. (2008 



RCO^ 
HCN 



5 (12) 
2 (13) 



One reason is that our disk is not meant to be specific 
to TW Hya even though we adopt the X-ray and UV 
spectrum of this source in our model. Another reason is 
that our column densities involve integration over the full 
vertical extent of the disk from the lower disk surface to 
the upper disk surface. In reality, the disk is likely to be 
optically thick, for 12 CO and H 12 CO + emission partic- 
ularly, and so the observed line emission is only tracing 
a fraction of the total column density. A third reason is 
that approximations for the disk chemical and physical 
structure are used for generating observed column densi- 
ties rather than a detailed model which allows for vertical 
and radial structure in the disk density and temperature 
and resulting molecular abundances. 

As stated previously, in order to determine how accu- 
rate our model is, one must perform full radiative trans- 
fer and compare our modelled line profiles with those 
observed. 

4. SUMMARY 

We have investigated the effects of photochemistry and 
X-ray ionization on the chemical structure of a proto- 
planctary disk surrounding a typical T Tauri star. We 
used a high-resolution complex physical model which 
takes into account irradiation by both the central star 
and the interstellar medium. We compiled a compre- 
hensive chemical network including a large gas-phase re- 
action set extracted from the UMIST Database for As- 
trochemistry (Rate06), gas-grain interactions (accretion, 
thermal desorption, cosmic-ray-induced desorption and 
photodesorption) and a grain-surface network. 

In previous work, we presented results from a model 
in which we approximated the photorates in the disk by 
scaling those in Rate06 (which assume the interstellar 
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radiation field) by the wavelength-integrated UV flux at 
each point. The X-ray ionization rate everywhere was 
also approximated using a power law to describe the 
energy-dependent cross section for ionization. Here, we 
recalculated both the photorates and the X-ray ioniza- 
tion rate, in the first case taking into account the unique 
UV wavelength spectrum at each point, and in the sec- 
ond, the unique X-ray energy spectrum and the explicit 
elemental composition of our gas. We also added the 
direct X-ray ionization of elements. 

We found that the recalculation of the photochemistry 
has a much larger effect on the disk chemical structure 
than that for the X-ray ionization. Concentrating on 
those molecules which have been observed or searched 
for in disks at either (sub)mm or infrared wavelengths, 
we find that the species most sensitive to the photochem- 
istry are OH, HCO+, N 2 H+, H 2 0, C0 2 and methanol. 
We also find the radicals, CN and C 2 H, are also af- 
fected although not to the same extent as those listed 
above. The recalculation of the photochemistry affects 
each molecule in a different manner indicating the in- 
herent non-linearity of the chemistry. Molecules affected 
throughout the disk extent include HCO + , N 2 H + and 
OH i.e., predominantly molecular ions and radicals. Wa- 
ter and methanol are mainly affected beyond a radius of 
around 1 AU, whereas C0 2 is altered beyond approxi- 
mately 10 AU. The differences in the behaviour of the 
saturated molecules can be attributed to the location of 
the molecular layer in each case, with methanol and wa- 
ter residing in a layer higher in the disk than C0 2 which 
in turn is linked to the desorption temperature of both 
molecules, with water more strongly bound to the grain 
surfaces than carbon dioxide (see Table [2]). 

In general, we find the change in the photochemistry 
leads to a depletion of molecules in the inner disk (within 
1 AU) and an enhancement in the outer disk (beyond 
1 AU). We conclude that, in the outer disk, our orig- 
inal approximation was overestimating the photorates, 
whereas, in the inner disk, where the stellar irradiation is 
strongest, we were underestimating the photorates. We 
should note that the photochemistry also indirectly af- 
fects the abundances of molecules via gas-phase chem- 
istry. In models where we recalculate the photorates, we 
find differences in the location of the H/H 2 and C + /C 
transition regions with both ratios generally decreasing 
in the molecular layer in models where we recalculate the 
photochemistry. Neutral molecules such as water and 
C0 2 are effectively destroyed by reaction with C + and 
in regions of high temperature by reaction with H atoms. 
Also, the molecular ions, HCO + and N 2 H + are efficiently 
destroyed via dissociative electron recombination. 

Concerning the X-ray ionization rate, we find the only 
molecule especially affected is N 2 H + . Residing in a differ- 
ent region of the disk to those molecules discussed thus 
far, we see the effect of X-ray ionization in the outer 
disk, in particular, where the abundance and distribu- 
tion of N 2 H + is decreased. Since the chemistry of N 2 H + 
is linked closely to the ionization of molecular hydrogen 



in the disk, this is unsurprising, however, one conclusion 
we can draw is that observations of N 2 H + emission from 
disks in conjunction with observations of those molecules 
unaffected by X-rays e.g., CO, could give us information 
on the effects of X-ray ionization in the outer regions of 
protoplanetary disks around classical T Tauri stars. 

The ionization fraction in the disk is affected only 
marginally by the recalculation of the photochemical and 
X-ray ionization rates. We find photochemistry affects 
the electron abundance in the disk surface whereas X- 
rays affect the disk midplane and outer disk region, since 
X-rays are able to penetrate the disk more effectively 
than UV photons. 

A calculation of the magnetic Reynolds number and 
ambipolar diffusion parameter, everywhere in our disk 
allowed us to determine the location of a possible sig- 
nificant dead zone in the disk midplane extending to 
« 200 AU where accretion will likely be suppressed. This 
has implications on both the physical and chemical struc- 
ture since a significant source of heating in the disk mid- 
plane close to the star is viscous dissipation due to ac- 
cretion flow. Our calculations suggest the accretion rate 
in the outer regions of our disk model is only » 70 % 
that used to determine the disk physical structure. In 
the inner disk, this value falls to ~ 1 %, however, we 
have neglected alternative sources of ionization which 
may dominate in the inner disk e.g., the thermal ion- 
ization of alkaki models, which should be considered in 
future models. 

We conclude that a wavelength-dependent treatment 
of the photochemistry in protoplanetary disks is neces- 
sary since it significantly affects the chemical structure. 
We also determine that existing approximations of the 
X-ray ionization rate in such objects are sufficient, af- 
fecting only the abundance and distribution of N 2 H + . In 
order to directly compare our model results with obser- 
vation and to determine the effects of the various chemi- 
cal processes investigated in our work on observable line 
emission we have computed the radiative transfer and 
modelled molecular line emission from the disk both at 
the resolution of existing facilities and at the expected 
resolution of ALMA. We report the results from these 
calculations in a subsequent publication (C. Walsh et al. 
2011, in preparation). 
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